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1. Introduction 

Recently, observational data of thermally emitting isolated neutron stars (NSs) confirm that 
most of them have magnetic fields larger than 10 13 G. Therefore, a reliable treatment of the thermal 
evolution must not avoid the inclusion of the effects produced by the presence of high magnetic 
fields. 

The non-uniform distribution of the surface temperature of isolated NSs seems to be con- 
firmed by the analysis of observational data (see reviews |[T]] and [^]). The mismatch between 
the extrapolation to low energy of the fits to X-ray spectra, and the observed Rayleigh Jeans 
tail in the optical band {optical excess flux), cannot be addressed with an uniform temperature. 
Several simultaneous fits to multiwavelength spectra of RX J1856.5-3754 [g], RBS 1223 [g], and 
RX J0720.4— 3125 [Q] are explained by a small hot emitting area ~ 10-20 km 2 and an extended 
cooler component. 

It has been proposed that the non-uniform surface temperature distribution may be produced 
by crustal confined magnetic fields [||, ^]. Magnetic fields larger than 10 13 G limit the movement 
of the electrons (heat carriers) in the direction perpendicular to the field with the result that the 
thermal conductivity is highly suppressed, while remains almost unaffected along the field lines. 

For such large fields, the field decay through Ohmic dissipation and Hall drift processes is very 
efficient and the heat released in the crust (Joule heating) must be taken into account in the thermal 
evolution of a neutron star [§]. In this article we focused on the effects of field decay and Joule 
heating on the neutron star cooling. In particular, we compare our simulations with observational 
data of a sample of isolated NSs that are highly magnetized. 



2. Cooling of neutron stars with magnetic fields 

We have performed two-dimensional simulations by solving the energy balance equation that 
describes the thermal evolution of a neutron star (NS) 

dT 

C v ^-V-(k-VT) = -Q v + Qj, (2.1) 

where C v is the specific heat per unit volume, Q v are energy losses by v-emission, Qj the 
heat released by Joule heating, and it is the thermal conductivity tensor, in general anisotropic in 
presence of a magnetic field. In this equation we have omitted relativistic factors for simplicity. 
A detailed description of the formalism, the code, and results can be found in [^]. The geometry 
of the magnetic field is fixed during the evolution. As a phenomenological description of the field 
decay, we have assumed the following law 

d — d exp(-f/T 0hm ) 

1 + (Tohm/ T Hall)(l ~ exp (-f/Tbhrn)) 

where B is the magnetic field at the pole, Bo its initial value, Tohm is the Ohmic characteristic time, 
and Tnaii the typical timescale of the fast, initial Hall stage. In the early evolution, when t <C Tbhm> 
we have B ~ Bq{\ +?/THaii) 1 while for late stages, when t > Tohm> B ~ Bo&xp(— f/Tohm)- This 
simple law reproduces qualitatively the results from more complex simulations |Q] and facilitates 
the implementation of field decay in the cooling of NSs for different Ohmic and Hall timescales. 
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3. Joule heating governing the cooling 

Given a 1.35 M Q neutron star model with a crustal confined magnetic field as in we have 
varied the parameters that describe the typical timescales for Ohmic dissipation and a fast initial 
decay induced by the Hall drift. 





Figure 1: Cooling of strongly magnetized NSs with Joule heating for Bo = 5 x 10 G. Upper panel shows 
T[, vs. t at the pole (left) and at the equator (right). Lower panel shows T s normalized to its value at the pole 
vs. 6. Three pairs of Joule parameters (Tbhm, Tkaii) are snown: (10 6 yr, 10 3 yr) with solid lines, (10 6 yr, 10 4 yr) 
with dashed lines, and (10 7 yr, 10 5 yr) with dotted dashed lines, respectively 

In Fig. [I] we show the cooling curves for different values of Joule parameters (Tbhm, ^Maii): 
(10 6 yr, 10 3 yr), (10 6 yr, 10 4 yr), and (10 7 yr, 10 5 yr) represented by solid lines, dashed lines, and dash- 
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dotted lines, respectively. For comparison, the thin dashed lines show the evolution with constant 
field for the same initial field B = 5 x 10 14 G. 

It is first to notice that there is a large effect of the field decay on the temperature at the bottom 
of the envelope 7/,: as a consequence of the heat released, it remains much higher than in the case of 
non-decaying magnetic field. The strong influence of the field decay is evident for all parameters 
chosen. The temperature of the initial plateau is higher for shorter THaii, but the duration of this 
stage with nearly constant temperature is also shorter. When t = THaii, B has decayed to about 
1 /2Bq and 3/4 of the initial magnetic energy has been dissipated. After t = Tnaii, Tb drops due to 
the transition from the fast Hall decay to the slower Ohmic decay. 

The insulating effect of tangential magnetic fields is twofold. First, in the absence of additional 
heating sources, it decouples low latitude regions from the hotter core resulting in lower tempera- 
tures at the base of the envelope. Second; if there is heat released in the crust, it prevents the extra 
heat to flow into the inner crust or the core where it is more easily lost in the form of neutrinos. 
Our simulations with Joule heating show systematically a hot equatorial belt at the crust-envelope 
interface. However, as discussed in [§], the inverted temperature distribution at the level of the 
crust is not necessarily visible in the surface temperature distribution because it is filtered by the 
magnetized envelope. An analysis of the angular temperature distribution given in the lower panel 
of Fig. [j] shows the development of a middle latitude region hotter than the pole at relatively late 
stages in the evolution (t ~ 10 4 , 10 5 yr). This hotter area is found with a wide range of parameters, 
and it would have implications on the light curves of rotating NSs, that will differ substantially 
from the light curves obtained with a typical hot polar cap model. 



4. Comments on the spin-down age of NSs 



Another aspect that should be reconsidered when we try to fit cooling curves to observations 
is that for many objects the age is calculated from the measurements of the rotational period P 
and its derivative P. The spin down age, t s d = P/2P, is derived assuming that the lose of angular 
momentum is entirely due to dipolar radiation from a constant (in time) magnetic dipole. In the 
case of a decaying magnetic field, ? sc j can seriously overestimate the true age t. A simple algebra 
shows that in the case of purely Ohmic decay 



2 V T Ohm/ 



(4.1) 



In the case that a Hall-induced fast decay also occurs, Eq. (2.2) results in a large correction of ? S( j 
as follows: 



hd = THall f(t) e 2t '^ 



m 



-'/ T Ohir 



Tbhm 



m^f{t) 



(4.2) 



where /(?) = 1 + — g-'/ T ohm) This relation gives t s d S> t by several orders of magnitude for 

t S> THaii, as shown in Fig. ^. Therefore, the cooling evolution time should be corrected according 
to the prescription for the magnetic field decay in order to compare with the observations properly. 



A detailed comparison with observational sources is presented in [10] and is summarized in the 
next section. 
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t(yr) 



Figure 2: Spin down age (f s( j) vs. true age (t) for a model with magnetic field decay. Results for two rates 
(fohm, iHaii) arc shown: (10 6 yr, 10 3 yr) with solid lines and (10 6 yr, 10 4 yr) with dashed lines. The grey solid 
line represents f s( j = t, and the short dashed line the purely Ohmic decay with Tbhm = 10 6 Y r - 



5. Comparison with observations 

We compare in Fig. |3| our simulations with observational data of NSs covering about three 
orders of magnitude in magnetic field strength: from radio-pulsars (B ~ 10 12 G) and isolated radio- 
quiet NSs (B ~ 10 13 G) to recent magnetar candidates (B ~ 10 14 ~ 15 G). The sources considered 



here are listed in Table 1 of [10], with the corresponding references. 

For most NSs, B is estimated by assuming that the lose of angular momentum is entirely due to 
dipolar radiation. The dipolar component is Bd = 3.2 x iO 19 (PP) 1 / 2 G, where P is the spin period 
in seconds, and P is its time derivative. In order to work with an homogeneous sample, we have 
included in the comparison only those objects for which P is available and the quoted magnetic 
field is B c j and discarded those sources for which B is inferred by other methods 1 . 

The reported temperatures are in most cases blackbody temperatures, except for low field 
radio-pulsars for which we take the temperature consistent with Hydrogen atmospheres following 



the criteria in |JL lfl. Nevertheless, there are some objects in which the estimate is an upper limit for 
the thermal component, like the Crab pulsar. This is also the case for some magnetars, which show 
large variations in the flux in the soft x-ray band on a timescale of a few years, indicating that the 
thermal component must be measured during quiescence and that the luminosity during their active 
periods is a result of magnetospheric activity. 

The age of a NS is subject to a large uncertainty, but it can be estimated by the spin-down age 



'For a few radio-quiet isolated NSs B can also be estimated assuming that observed x-ray absorption features are 
due to proton cyclotron lines, but this gives the surface field, which is usually larger than the external dipolar component 
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Figure 3: Cooling curves with corrected spin down age. Observational sources correspond to Table 1 of [[!(]] . 
Symbols identify sources with the same order of magnetic field: squares for magnetars (AXPs and SGRs, 
with B ~ 10 14 ~ 15 G), triangles for radio-quiet isolated NSs (with B ~ 10 13 G), and circles for radio-pulsars 
(with B ~ 10 12 G). Open circles denote temperatures obtained from fits to Hydrogen atmospheres. 



(tsd = P/2P), provided that the birth spin rate far exceeds the present spin rate and B is considered 
constant 2 . As shown in Sec. ||, if one considers magnetic field decay, t s d seriously overestimates 
the age t of the simulations. Therefore, we transform the cooling curves to plot T s vs. t st j, to include 
the temporal variation of the magnetic field. 

In Fig. |5| our results show that the effect of high magnetic field (Bo — 10 14 ~ 15 G, solid lines) 
in the cooling is important from the very beginning of the NS evolution. The temperature reached 
is increased up to a factor of 5 in comparison with a non-magnetized model and can be kept nearly 
constant for about 10 4 years. The effect of Joule heating is very significant and may help to explain 



why magnetars are so hot [12]: the high temperatures in the early epoch result in higher electrical 
resistivity and in an faster magnetic field dissipation that releases the heat in the crust. In this 
picture, the thermal evolution of radio-quiet, isolated NSs could be represented either by NSs born 
with intermediate fields in the range of Bo = 10 13 -10 14 G (dashed lines) or by magnetars in which 
the field has already decayed in a timescale of sa 10 5 ~ 6 years. For intermediate field strengths, 
the initial effect is not so pronounced but the star can be kept much hotter than non-magnetized 



"For some cases an independent kinematic age is available, which does not necessarily coincide with t sc [. 
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NSs from 10 4 yr to 10 6 yrs. For weakly magnetized NSs, radiopulsars with B ~ 10 12 G, the effect 
of the magnetic field is small (dashed dotted lines) and they can satisfactorily be explained by 
non-magnetized models, with the exception of very old objects (t > 10 6 years), as discussed in 

B 



6. Conclusions 

From the results presented here we conclude that the thermal evolution of a magnetized NSs 
is strongly affected by the presence of magnetic fields. Therefore, studies aimed to disentangle 
properties of NSs interior (e.g. EoS, neutrino processes, etc.) for objects with B > 10 13 , through 
cooling curves should not neglect the (dominant) magnetic field effects. A first step towards a cou- 
pled magneto-thermal evolution has been given in this work and future investigations will consider 
a consistent evolution including the evolution of the magnetic field geometry. 
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